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Abstract In this paper we study spectral function for a nonsymmetric dif¬ 
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considered, and for each case we prove by Marchenko’s method that, to the 
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1 Introduction 

As a very essential mathematical problem, Weyl-Stone eigenfunction ex¬ 
pansion |291 132] in which the key role is spectral function for singular self- 
adjoint second order linear differential operators has been studied deeply by 
many renowned mathematicians such as K. Kodaira m, N. Levinson dZI> B. 
M. Levitan [18], E. C. Titchmarsh [30] and K. Yosida [34] and so on. This the¬ 
ory for the singular operators may be derived as a limiting case of the classical 
Sturm-Liouville expansion theorem for the regular operators, where the Parse- 
val equality for the regular operators plays a very important role for the proofs. 
Similar ideas can also be applied to singular self-adjoint first order systems, for 
example, the Dirac operators dans]. For general theory of eigenfunction ex¬ 
pansion for self-adjoint and regular non-self-adjoint operators in Hilbert space, 
we refer to mmm- For multidimensional cases, see, e.g., m ■ Moreover, 
a two-fold spectral expansion in terms of principal functions of a Schrodinger 
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operator has been derived in pQ. Recently, K. Kirsten and P. Loya m have 
obtained some interesting results on the spectral zeta function for a Schrodinger 
operator on the half line. 

However, to the author’s knowledge, for singular nonsymmetric differential 
operators, there are few results on eigenfunction expansion. The limiting ap¬ 
proach for self-adjoint case can not be applied even for very simple case of 
nonsymmetric differential operators, since in general the corresponding regular 
spectrum has irregular behavior on the complex plane. In order to extend expan¬ 
sion theory to general case, V. A. Marchenko mm established an excellent 
method in dealing with the singular Sturm-Liouville operator with complex- 
valued potential. In this paper, inspired by the idea of V. A. Marchenko, we are 
going to establish expansion theorem in two cases for a singular nonsymmetric 
differential operator, where the key is to prove the existence of the correspond¬ 
ing spectral function. Our results can be extended to 2 n x 2 n systems, and 
for simplicity we here will only consider the case of n = 1. For the regular 
case of this nonsymmetric differential operator, recently we have obtained some 
results on inverse spectral problems with applications to inverse problems for 
one-dimensional hyperbolic systems, see It is well known that for many 

differential operators there are intrinsic relations between their spectral func¬ 
tions and the corresponding Weyl functions (often called m- functions), and for 
the recent interesting results on Weyl functions see, e.g., 

For the asymptotic behavior of spectral functions for elliptic operators we refer 

to mmm 

In this paper we consider boundary value problems generated by a nonsym¬ 
metric differential operator on the half line 0 < x < oo: 

(A P ip)(x) := B^-(x) + P(x)ip(x) = \ip(x) 


where B = 


both matrix-valued function P = 


P li 
P 21 


P 12 
P22 


(C 1 [0,oo)) 4 and parameter A are complex-valued. It is directly checked that 
the adjoint operator of Ap in some suitable Hilbert space is — Bf^ + P T (x) and 
consequently Ap is nonsymmetric. Here and henceforth, c denotes the complex 
conjugate of c and - T denotes the transpose of a vector or matrix under con¬ 
sideration. Here we point out that the spectrum problem for Ap with compact 
matrix-valued function P has been studied in Ell¬ 
in order to describe our results properly, we first give some information on 
distributions and we refer to j2Tj for more details. Let K 2 (0,oo) denote the set 
of all square integrable functions in (0,oo) with compact support. For a > 0, 
we set K 2 (0,oo) = {/ € K 2 (0, oo) : f(x ) = 0 for x > a}. The entire function 
e(p) is called the function of exponential type if |e(p)| < C'exp(cr|Imp|) where 
the positive constants C and a depend on e(p). Moreover, the index 


a P = lim r 


-l 


In ( max |e(p)| 


\|p|=r 
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is called the type of entire function e(p). Let linear topological space Z be 
the set of all entire exponential type functions integrable on the real line. The 
sequence e n converges to e in Z if limn^oo \e n (p) — e(p)|dp = 0 and the 
types cr„ of the functions e n (p) are bounded: sup a n < oo. The set of all linear 
continuous functionals defined on the test space Z will be denoted by Z' whose 
components are called distributions (generalized functions). The sequence D n 
converges to D in Z' if lim^oo < D n , e(p) >=< D, e(p) > for all test functions 
e G Z. 

In this paper we consider two cases of the coefficient matrix P. The first case 
is special and will be described as follows. Let P be a continuously differentiable 
matrix-valued function satisfying BP = PB and p be a complex constant. 

Here it is easy to see that P is of the form ^ ^ ) Consider the following 

boundary value problem 


B^-(x) + P(x)p(x) = Xp(x), 0 < x < oo, 
ax 


m = 


cosh p sinh p 
sinh p cosh p 


( 1 . 1 ) 


Let ip = < p(x, A) be the solution to (11.111 and 

¥>[i] = I and ¥>[ 2 ] = 


( 2 ) 
V *>[1] 


rm 


be the first and the second column vector of the matrix <p, i.e., p = {p[\\ P[ 2 \)- 
Similarly we denote the matrix inverse of p by = p^ 1 = (i/>m ipn])- Now for 

/ = ( f f w ) G (^ 2 (0,oo)) 2 , g = ( ^ ) G (Z/ 2 (0,oo)) 2 

where (L 2 (0, oo)) 2 denotes the product space of L 2 { 0, oo), we set 

POO POO 

U}k f(p)= f T (x)i> [k] (x,ip) dx, T] k g {p) = / pf k] (x, ip)g(x)dx {k = 1,2), 

J 0 Jo 

where i = \/—l, p £ R. Then we have the first main result of this paper. 
Theorem 1. It holds for the boundary value problem m that 

poo 1 2 fOO 

/ f T (x)lKx) dx=—J 2 u k f(pH(p)dp. (1.2) 

JO Z7r k=1 J-oo 

Moreover, for f G (K 2 (0,oo)) 2 with uj k (p),r] k (p) £ Z (k = 1,2), the following 
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expansion formula holds: 


f(x)= 


1 

2n 



W/(pV[fc] ( X,ip)dp 


1 

27r 



V f(p)^[k] OM/?)dp. 


(1.3) 


We often call (11.21) (or (11.71) 1 the Marchenko-Parseval equality which means 
that a spectral function exists in corresponding boundary value problem. His¬ 
torically, the concept of spectral function came from the classical theory of Weyl. 
Theorem 1 implies that is a spectral function corresponding to problem 
EH) with P satisfying BP = PB , which is the same as the case of P = 0. Here 
and henceforth E denotes the 2x2 unit matrix. 

For general matrix function P £ (C' 1 [0,oo)) without the constraint BP = 
PB , we also can show the existence of the corresponding spectral function. 
More precisely, let Q be a 2 x 2 matrix satisfying QB + BQ = B and Q 2 = 
Q. It is seen by simple computation that there exists matrix Q satisfying the 
above conditions, and the simplest one is Q =diag(l,0). It follows easily from 
detH = 1 that detQ = 0. Consider the following boundary value problems 


! B-^-(x) + P(x)ip(x) = Xtp(x), 0 < x < oo, 
ax 

¥>(0) = Q, 


(1.4) 


and 


— ^-{x)B + ip[x)P[x) = \ip(x), 0 < x < oo, 
dx 

£(0) = Q- 


(1.5) 


Denote the solutions to problems EH and EH by <p(x, A) and ip(: r, A), respec¬ 
tively. For all 2x2 matrices f,g £ (i 2 (0,oo)) , we set 



f{x)<p(x,ip)dx, 



p(x,ip)g{x)dx, 


( 1 . 6 ) 


where i = y/—l, p € R. Then we have another main result of this paper. 


Theorem 2. To the problems Q and EH there corresponds a distribution¬ 
valued spectral function D = (Dki) 1<k t<2 such that D = QDQ, D k i £ Z' and 


f(x)g{x)dx = / $f(p)D(p)$ g {p)dp. 


(1.7) 
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Moreover, for f £ (K 2 (0,oo)) 4 with $/(p),<f>y(p) £ Z 4 , the following expansion 
formula holds: 


/ oo /*oo 

®f(p) D (p)<P(x>ip)dp= / ‘P(x,ip)D(p)$f(p)dp. (1.8) 

-oo J —oo 


Although Theorem 1 and Theorem 2 have shown the existence of spectral 


function for the singular nonsymmetric differential operator in two cases, we 
here point out that the uniqueness of spectral function for the operator does 
not hold generally, which is the same as that for Sturm-Liouville operators (see, 
e.g., |19]1. Moreover, since the spectral function is distribution-valued, it is not 
a measure in general, which is different from the case of self-adjoint Sturm- 
Liouville operators. Besides, given singular nonsymmetric differential operators 
with general P, it is still an open problem to prove the existence of spectral 
functions under general boundary conditions. On the other hand, it is inter¬ 
esting to investigate the corresponding inverse problems, namely, given spectral 
functions or Weyl functions, find the differential operators. See [5] for the clas¬ 
sical inverse problem to determine the potential of the Sturm-Liouville operator 
from its spectral function and [Ij for determination of singular differential pen¬ 
cils from the Weyl function. Theorem 1 has implied that the uniqueness does 
not hold generally for the inverse problems, and we need impose other assump¬ 
tions for uniqueness. In a forthcoming paper we will study the inverse problems 
for the singular nonsymmetric differential operator. 

The paper is composed of four sections. In Section 2 we establish transfor¬ 
mation formulae for our boundary value problems. Sections 3 and 4 are devoted 
to prove Theorem 1 and Theorem 2 by transformation formulae, respectively. 


2 Transformation formulae 



( 2 . 1 ) 



( 2 . 2 ) 


and 



(2.3) 


Moreover let us put 



(2.4) 
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Here we remark that P(Pl,P 2 )(0) = E, P(Pi,P 2 )(x) = P _1 (P 2 , Pl)(x) and 
R Pfi, —P 2 ^J (x) = R (P 2 ,Pi) (x). Let M 2 (C) be the set of all 2x2 complex¬ 
valued matrices. We first prove the following lemma. 


Lemma 2.1. For any A £ C, Q € M 2 (C) with det Q = 0 and Pj £ (C 1 [0, oo)) 4 
{j = 1,2), let ipj = ipj(x, A) satisfy 


+ Pj(x)(fj(x) = \tpj(x), 0 < x < oo, 
<Pj( 0 ) = Q- 


(2.5) 


Then there exists a unique K(Pi,P 2 ;Q) = (Kki(Pi, P 2 ] Q)) 1<fc j <2 £ (C 1 ^)) 4 
independent of A such that for 0 < x < oo and all A £ C 

T 2 (x, \) = R(Pi, P 2 )(x)tpi(x, \) + f K(P 1 ,P 2l Q)(x,y)(p 1 {y,X)dy. (2.6) 

Jo 


(transformation formula) 

Here R(P±, P 2 )(x ) is defined by 

Moreover, the kernel K(P\, P 2 ;Q) is the unique solution to the following 
problem of first order system 


+P 2 (x)K(P 1 ,P 2 ,Q)(x,y) - K(P 1 ,P 2 ;Q)(x,y)P 1 (y) =0, ( x,y ) S H. 

(2.7) 

I\(Pi,P 2 ; Q)(x, 0)BQ = 0 (0 < x < oo). (2.8) 


K(P 1 ,P 2 -Q)(x,x)B- BK(P 1 , P 2 ; Q) (x, x) 

= B dR ^ Pl ' P2) (x) + P 2 (x)P(Pl, P 2 )(x) - P(Pi,P 2 )(x)Pi(x) 


dx 


(2.9) 


(0 < x < oo). 


Proof. We prove this lemma by the idea used in [33]. Since Pj £ (C^fO, oo)) 4 

(j = 1,2), it can be verified directly that, if K(Pi,P 2 ,Q) £ (C 1 (H)) 4 is the 
unique solution to problem (12.71) ^ (12.91) . then (12.61) holds. Therefore, it is suf¬ 
ficient to prove the existence and the uniqueness of the solution to problem 
(12.71) ^ (12.91) for each Pi,P 2 S (C ,1 [0,oo)) 4 . 

For clarity, we reduce the proof to a special case. By the condition det Q = 0, 
we may assume that a complex constant c exists such that q 2 = cqi where q \, q 2 
are the first column vector and the second one of Q , respectively. Then it is 
sufficient to prove the existence and the uniqueness of the solution to problem 
(12.71) ^ (12.91) in the case <Pj(d, A) = q\, since problem (12.51) is linear. Moreover, 

since a complex constant c* exists such that q± = c* ^ g inhp ) w ^ iere M G C, 
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cosh fi 
sinh/i 

the the solution to problem (I2.71> ~ (|2.9[1 by K(Pi,P 2 ,y)(x,y), and (12.811 has the 
following form: 


J . In this case, we denote 


it can be reduced to the case ipj (0, A) = 


f K 12 (P 1 ,P 2 ,y.)(x, 0) = -tanh/z K 11 (P 1 ,P 2 ,fj,)(x,0), 
\ K 22 (P 1 ,P 2 ,y)(x,0 ) = -tanh/z K 21 (P 1 ,P 2 ,/j)(x,0). 

If we set 

( L\(x, y) = K 12 (P 1 ,P 2 ,y)(x,y) - K 21 (P 1 ,P 2 ,y)(x,y), 
I L 2 (x,y) = K 11 (P 1 ,P 2 ,y)(x,y) - K 22 (P 1 , P 2 , y,)(x,y), 
I L 3 (x,y) = K 11 (P 1 , P 2 , y)(x,y) + K 22 (P 1 , P 2 , y,)(x,y), 
{ L 4 (x,y) = K 12 (P 1 ,P 2 ,y)(x,y) + K 21 (P 1 ,P 2 ,y)(x,y) 


( 2 . 10 ) 


( 2 . 11 ) 


and L = L(x,y) = (Li(x,y),L 2 (x,y),L 3 (x,y),L 4 (x,y)), then we can rewrite 


(12.71)^(12.91) as follows: 




dL k (x,y) dL k {x,y) 

ax a y =fk{x ' y ’ L) 

((x,y) e fi, 

A =1,2), 

(2.12) 

9L k (x, y) dL k {x,y) 

ax + a y =*<*■»• L > 

(0*5 2/) G fi, 

k = 3,4), 

(2.13) 

L k (x,x) = r k {x) 

(0 < x < oo, 

fc = l,2), 

(2.14) 

f L 3 (x,0 ) = sinh(2/z)Li(a:, 0) + cosh(2/u)L 2 (x, 0) 

[ L 4 (x,Q) = — cosh(2^)Li(a;, 0) — smh{2fj)L 2 {x, 0) 

(0 < x < oo), 

(2.15) 


where f k (x,y,L) = ^J2m=i( a km(y) + b k m(x))L m (x,y) (1 < k < 4), here 
a>km{y) , b km (x ) (1 < fc,m < 4) are linear combinations of two elements of 
the matrix functions P\(y) and P 2 (x) respectively, and r k G C' 1 [0,oo) (k = 1,2) 
are dependent only on P\ and P 2 . 

Integrating (12.128 . (12.131) with (12.141) and (12.151) along the characteristics 
x + y = const, and x — y = const, respectively, we obtain the following integral 
equations: 

Lk{x,y)= f k {-s + x + y, s, L)ds + r k ( X ^ V ) {{x, y) G fi, k = 1,2), 

(2.16) 

and 

ry 

L k (x,y) = / f k (s + x - y,s,L)ds 

Jo 

x — y 

+ / {a k fi{~s + x - y, s, L) + Pkf 2 (s + x-y,s, L)} ds 
Jo 

,x — y. _ ,x — Vx 

+«fcri(^-^) + /3kr 2 (—2^) 


(0*5 2/) G fi, fc = 3,4) , 


(2.17) 
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where az = sinh(2/x), fa = cosh(2^t) and a .4 = — cosh(2 y), /? 4 = — sinh(2 y). 

The unique solution L £ (C 1 (n)) 4 to (12.161) and (12.171) can be obtained by 
the iteration method. In fact, setting 

L k\ x ,y) = 0 (( X ,V) e to, 1 < k < 4) , 

x±y_ 

l[ n \x,y) = J fk (-s + z + y,s,L (n_1) ) ds + r fe (^-^) 

({x,y) £ Cl,n > l,k = 1,2), 

and 

L i n) ( x ,y) 

= J fk (s + x-y,s, L (n_1) ) ds 



|a fc /i (-s + x- y, s,T ( " + j3 k f 2 (-s + x - y,s,L (n 4) j j ds 


+a k ri{ 


x-y. 


■P k r 2 (^) 


2 ' ' 2 

((z,y) e 0,n > 1, k = 3,4), 

we can obtain by induction the estimates for each n > 1 

C n ~\x) 


L^{x,y)-L^{x,y) <u{x) 


(n — 1)! 


((x,y) £Sl,l< fc<4), (2.18) 


where 


u/(x) = (| sinh(2/x)| + | cosh(2/Li)| + 1) max (|ri(s)| + |r 2 (s)|) 

0<s<x 


and 

l 2 

C(x) = (| sinh(2/x)| + | cosh(2^i)| + 1) x max - V (|P lifc/ (s)| + \P 2 . k i{s)\) ■ 

0<s<x Z z ' 

- k,l =1 


Thus L k (x,y) = lim„_ ! . 00 L^\x,y) (1 < k < 4) exist uniformly for (x,y) £ ft 
and we see that L k {x,y) (1 < k < 4) satisfy (12.161) and (12.171) with the bound 
\L k {x,y)\ < w(x)exp{cr(x)). 

Moreover, differentiating (12.161) and (12.171) with respect to x and y, we can 
similarly obtain by induction the following estimates 


dL { k \x,y) 

dx 


dL { k n l \x lV ) 
dx 


< Z(x) 


(n — 1)! 


((x,y) 6fl,l<K4), 

(2.19) 
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dL^\x,y) 

dy 


where 


dL ( k l) [x,y) 
dy 


< £0) 


(«-!)! 


(( x,y ) G O,1 < k < 4) , 

( 2 . 20 ) 


= - (| sinh(2/i)| + | cosh(2/i)| + 1) 

x I ™ a .? d r i( s )l + l r 2( s )l) + ^(i) exp(C(x)) 
k0<s<x Z 

2 

X o<s<x (l^ > i>fei( s )l + l-P2,fci( s )l + (IA , ,fc;( s )l + l-f2,fe;( s )l) x ) }• 

- s - x k,i =l 

Therefore, it follows from (12.191) and (12.201) that L G (C ,:L (fl)) 4 . The uniqueness 
of the solution to (12.7IWfl2.9l) is shown by (12.181) . □ 


Corollary 2.2. For j = 1,2, let ipj be the solution to problem 11.11) with 
P = Pj G (C 1 [0, oo)) 4 satisfying PjB = PPj. TTien 1/ie following transforma¬ 
tion formula holds: 

P 2 {x,X) = R(P 1 ,P 2 )(x)tpi(x,\) (2-21) 

where R(P\, P 2 )(x) is defined by 

Corollary 2.2 follows from the fact that K(Pi, P 2 , p) = 0, which can be de¬ 
rived easily by observing that the right hand side of (12.91) is 0 (in this case the 
condition det Q = 0 is not necessary). Or one may directly verify (12.211) . Here 
we omit the details. 


Corollary 2.3. 

Let S and S be the solutions corresponding to P = 0 in O and \1.5\) . respec¬ 
tively. Then the following transformation formulae hold. 

(1) For problem o we have 

S(x,ip) = R(P,0)(x)(p(x,ip) + f K(P,0;Q)(x,y)ip(y,ip)dy (2.22) 

Jo 

where the kernel K(P,0;Q) G (C 1 (H)) satisfies the equation 


BK X (P, 0; Q)(x, y) + K y (P, 0; Q)(x, y)B - K(P , 0; Q)(x, y)P(y ) 
as well as the conditions for 0 < x < 00 


0, (x,y) G n 
(2.23) 


K(P, 0; Q)(x, 0 )Q = K(P, 0; Q){x, 0), (2.24) 


and 


K(P , 0; Q){x, x)B-BK(P , 0; Q){x, x) = BR'(P, 0 ){x)-R{P, 0){x)P(x). (2.25) 
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(2) For problem U.5 \ ) we have 


S{x, ip) = ip{x , ip)R{ 0, P)(x) + [ <p{y, ip)K T (-P T , 0; Q T ){x , y)dy (2.26) 

J o 

where the kernel Jv T (—P T , 0; Q T )(x, y) satisfies 

QKT(-W,0-,QT)(x,0) = K T (-W,0-,CF)(x,0). (2.27) 


Proof. (1) is obvious, since detQ = 0 and then Lemma 2.1 can be applied. 
Here (12.241) follows from (12.81) . BQ = B — QB and B 2 = E. Now we prove (2). 
Note that by m the function <p(x, ip) statifies 

! B (x) — P T (x)p T (x) = ipip T (x), 0 < x < oo, 
dx 

Wio) = W- 


Then one will obtain (12.261) by (1) if he notices the following fact: f?(0, P){x) = 

R(-W,0)(x). □ 

Since the solutions to the boundary value problems with P = 0 are entire in 
A, it follows easily from the transformation formulae that 

Corollary 2.4. For each fixed x, all solutions to the boundary value problems 
under consideration are entire in A. 


3 Proof of Theorem 1 

We divide the proof of Theorem 1 into four steps as follows. 


First step. We first construct a regular spectral function. Let S denote the 
solution of m corresponding to P = 0. Set 

p = — iX and v = —ip. 


It is easy to see that 


S = S{x, A) 


cosh(Aa: + p) sinh(Ax + p) 
sinh(Ax + p) cosh(Ax + p) 


f cos (px + v) 
i sin(px + v) 


i sin(px + v) 
cos (px + v) 


(3.1) 
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and 


= S' -1 (a-, A) = 


(3.2) 


cos(px + v) —i sin(pcc + v) 

—i sin(px + v) cos (px + v) 

We choose two sufficiently smooth real-valued functions S n (x) and 7 <r(a;) subject 
to the following conditions: 


nOO 

/ S n (x)dx = 1, 
Jo 


/ 0 

S n (x) = 0 for x = 0 and x > — , 5 n ( x) > 0 for 0 < x < — , 

n n 

Jcr(x) = 1 for 0 < x < ( 7 , "/ a (x) = 0 for x > a + 1, (3-3) 

and it is obvious that S n (x) tends to the Dirac delta function S(x) as n —> oo. 
We set 


K(p) = (K, jm (p)) 




1 

27T , 


cos (px + v) —i sin(px + v) 
—i sin(pa; + v) cos(px + v) 


(3.4) 


xR(P,0)(x)6 n (x)E'y a (x) 


cos v i sin v 
i sin v cos v 


da’. 


Since the Fourier transform is a one-to-one mapping on the space of bounded 
continuous Lebesgue-integrable functions and R(P,0)(x)6 n (x)E'y a (x) is a con¬ 
tinuously differentiable matrix function with compact support, it is not hard to 
see that the matrix function D°{p) is bounded and Lebesgue-integrable on the 
real line — oo < p < oo. Hence the integral 


S(x,ip)D°{p) 


cos v —i smi/ 
— is\nv cosv 


dp 


converges absolutely. By Corollary 2.2 we have ip(x,ip)= R(0, P)(x)S(x, ip)= 
R~ 1 (P,0)(x)S(x,ip), which implies by the Fourier inverse transform that 


<p(x,ip)D°(p) 


cos v —* sin v 
—i sin v cos v 


dp = S n (x)E (0 <x<a). (3.5) 


Here and henceforth we repeatedly make use of the fact that two matrices P\ 
and P 2 in the form of ‘ 


are interchangeable: P 1 P 2 = PzPi- 


Second step. Next we will investigate the asymptotic behavior of the follow¬ 
ing matrix function as n —> 00 


UZ(x,y) = (yz tkl {x,y)j 


l<k,l< 2 


/ oo 

cp(x,ip)D^(p)ip~ 1 (y,ip)dp (0 <x,y< a). 

-OO 


(3.6) 
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It is easy to find that 


Un(x, o) = J <p(x, ip)D"(p) ^ 212 12 V ) dp = Sn ^ E (0 < X < a) 

and 


UZ(0,y) 


rt 

cosv 

i sin v \ 

Loo \ 

i sin v 

cosz^ J 


D n(p)v 1 (y,ip)dp 


(o <y< a). 


Now we show that U°(0,y) = 0 for all y > 0. Indeed, first one can see from 
(13.41) that 


K(p) 


n 

cos(px) 

—*sin(px) \ 

Jo V 

—i sin(px) 

cos (px) J 


x {i?n(P, 0 ){x)E + i?i 2 (P, 0)(x)P} s n (x)'y a (x)dx. 


Moreover, for any continuous scalar function u{x) with compact support and 
u(0) = 0, it follows easily from the theory of the Fourier cosine and sine trans¬ 
forms that 

( cos(px) —ism(px) \ , . / cos (py) —ism(py) \ . . 

( -(stater) cosO*) ) " {x) { -(state#) cos( ra ) ) 



= 0 . 

(3.7) 

Consequently, it follows from (13.71) and <p 1 (-,ip) = S 1 (-,ip)R(P, 0)(-) that 
Un(0,y)R(0, P){y) = o and hence U£(0,y) = 0, since R(0,P)(y) is invertible. 

On the other hand, by (13.61) it is easy to see that, for fixed n and a, U£(cr, •) 
is a bounded differentiable function on [0, a] and denoted by £„(•) for simplicity. 
Therefore, since by (O) we easily show that 

~~~t ^ B - ip~ l (x)P(x) = -ip<p _1 (x), 
ax 

the above argument implies that the functions 

/ N 

y(x,ip)D°(p)ip- 1 (y,ip)dp (3.8) 

-N 

are continuously differentiable and satisfy the equation 
dU dU 

B ~dx^ X,V ^ + ~d^^ x ' y ^ B+P ^ U ^ x ' y ^ ~ U ( x ’ y ) P ( y ) = 0 in ( 3 - 9 ) 

as well as the following conditions 

U(x, 0) = 5 n N(x)E (0 < x < a), (3.10) 
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U(0,y) = V nN (y), U(a,y) = E nN (y) (0 <y<a), (3.11) 

where IIo- = {( x,y ) el 2 : 0 < x,y < a}, the functions S n N,^nN and S n Ar sat¬ 
isfy the compatibility conditions and limjv^oo ^nN(x)= 5 n (x), limjv-x>o r n jv(y) = 
0 and lini/v->oQ S n jv(2/) = 3 n (y ). We should note that problem (13.91) . (13.101) and 
(13.111) can be rewritten as a symmetric hyperbolic system: 


dV . , 


0 E \ dV . 

— (x,y) + C{x,y)V(x,y ) = 0 m n CT , 


v P 0 

V(x,0) = J n jv(x)lt (0 < x < a), 

. V(0,y) = l* nN (y), V(a,y) =^ nN (y) (0 < y < a) 


where 


V(x,y) = 

( r n jv,ii(y) \ 


( U u (x,y) ^ 
Ui 2 (x,y) 


0 

U 2 i(x,y) 

5 12 — 

0 

\ U22{x,y) y 


l 1 / 


r„Ar(y) = 


iN 


(: V ) = 


r ra jv,i2(y) 

r„jv,2i(y) 

V r„jv,2 2 (y) / 

and C(x, y) is the following 4x4 matrix-valued function 

/ -Pi 2 (y) Pn(x) - P 22 (y) 0 


( H n jv,ii(y) \ 
S n jv,i2(y) 
3 n iV, 2l(2/) 

V 2 nN,22{y ) / 


Pn(z) - Pii(y) 
0 

P 21 W 


(3.12) 


Pi2(a:) 

-P 2 i(y) Pi 2 (z) 0 

P 2 l(®) -Pl2(y) P 22 (x) - P 2 2 {y) 

0 P 2 2 {x) — Pn(y) -P 2 i(y) 


Since BU^ N (x,y) = U% N ( x,y)B , a direct calculation shows that the symmetric 
hyperbolic system (13.121) is actually equivalent to the following normal hyper¬ 
bolic system 


— (x,y) = ^ q _ x J —(x,y) + c(x,y)v(x,y) in U a , 

v(x, 0) = 6 nN (x)h (0 < x < a), 

^2(0, y) = 2t>i(0, y) - r„jv,n(y) + 3r nA r,i2(y), 

„ v 2 (o",y) = 2vi(o-,y) - H n jv,n(2/)+33^,12(3/) (0 < y < a), 


(3.13) 
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where 


Un(x,y) - Ui 2 (x, y) 
Un{x,y) + U 12 (x,y) 


v{x, y ) 


and 

c(x,y) 


vi(x,y) \ 
V 2 (x,y) ) 




/ (P n - P 12 )(x) + (Pi 2 - P n )(y) (P 12 - Pu)(x) + (P 22 - P 2 l)(y) 
V (i’ll + P 12 M - (Pll + Pl 2 )(x) (P 22 + P 2 i)(y) - (P 11 + Pi 2 )(*) 


If we take the variable y as time, then it is not hard to verify that the classical 
Uniform Kreiss Condition holds, and hence from the well-known results of well- 
posedness on linear hyperbolic systems (cf. |5j and references therein) we see 
that ( 13 . 131 ) has a unique solution, that is, there exists a unique solution U° N (x, y) 
to problem ( 13 . 91 ) . ( 13 . 101 ) and ( 13 . 111 ) such that U° N (x,y) —¥ U%(x,y) as N —> 00 . 

On the other hand, if we set W% N (x, y) = U° N (x,y) — S n N(x — y)E for 
0 < x,y < a where 6 u n{x — y) = 0 for 0 < x < y < a, then W° N {x, y) satisfies 
the following equation 


B ~dx^ x,y ^ + ~d^^ x,y ^ B + p ( x ) w ( x ’y) ^ w ( x >y) p (.y) 

= S nN {x - y) ( P{y ) - P(x)) 


(3.14) 


and W[x,Q) = 0. It follows easily from the compatibility conditions that 
W„V0, y) 0 {N, n —> 00 ). Next wc will show 

W nN(^y)-^° (IV ,n —> 00 ). (3.15) 

In fact, it follows from ED, ED, ED and the transformation formulae 
( p(-,ip)= R(0,P)(-)S(-,ip)= R~ 1 (P,0)(-)S(-, ip), tp~ 1 (-,ip) = S~ 1 {-,ip)R(P, 0)(-) 
that 

s «(y) 


= P(0, P)(a)R(P, 0 )(y)R(P, 0)(a - y) 


X 



y,ip)K(p) 


costs — isiius 
—isinis cos is 


dp 


= R( 0, P)(<j)R(P, 0 )(y)R(P, 0)(cr - y)S n (a - y)E (0 < z < a). 
Thus, we have 


S„(y) - 5 n (a - y)E 

= 5 n (a - y)R{ 0, P)(a)[R(P, 0 )(y)R(P, 0 )(a - y) - R(P , 0)(a)] 


(3.16) 
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whence (13.151) follows easily. Consequently, by the wcll-posedness of symmetric 
hyperbolic linear differential equations, we have 

W° N { x i y) ~> 0 as N, n — > oo 

since S n n(x — y) —> 5(x — y) as N, n —> oo and hence the right hand side of 
(13.141) tends to 0. Therefore, for 0 < x,y < a 

UZ(x,y)^ 6 (x-y)E (noo). (3.17) 

Remark. There is another and simpler way to prove (13.171) in which it is not 
needed to consider m . The key idea is based on considering m, (ED and 
(13.161) with replacing a by x. We leave the details to the reader. 

Third step. We prove the Marchenko-Parseval equality m- Assuming that 
f,9 e (K2(0 , oo)) have compact support, we have by changing the order of 
integration that 

/-•oo /*oo 2 _ 

U n,kii x ^y)^ k) {x)f {l) {y)AxAy 

J 0 Jo 1U 7_1 


poo poo 2 / 

/ / E 

k,i =l \” 


I D n,jm(p) i p\j]{x,ip)i’[l! l] (y,ip)dp 

~°° j,m=l 


xg ( - k ){x)f ( - l \y)dxdy 


= df> D n,jm(P) 


j,m=l ' 


7*00 p OO 


f {l) (2/)^[m] (y> v)fi ,(fc) {x)dxdy 


/o Jo 


kd =1 


2 poo 

A ^ — J —OO 


j,m= 1 


Therefore, in view of (13.171) . we obtain by letting n —> oo that for any f,g £ 
(K2(0,oo)) 2 

p oo 2 7*00 

/ f T (x)g{x)dx = lim^ ^ / D* tjm {p)wY(p)r£(p)dp. (3.18) 
By the definition of (p) (see (13.41) ), we easily see that 


lim r%(p) = J-R(P, 0)(0) = J-E. 

n yoo Z7T Z7 T 


(3.19) 
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On the other hand, since the Fourier transform is a continuous mapping of 
L 2 (R) into L 2 (M.), it follows easily from Corollary 2.2 and the zero extensions of 
/ and g on R that both and r/^(p) belong to L 2 (R). Therefore, Combining 

(13.18j) and (13.191) and letting er —► oo, we can assert (11.21) by the boundedness 
of D°(-), the dominated convergence theorem and the fact that (K 2 (0, oo)) 2 is 
dense in (L 2 (0, oo)) 2 . 


Forth step. We prove the expansion m- First we assume that / £ (Co[0, oo)) 2 , 
where (Co[0,oo)) 2 denotes the product space of the set of all continuous func¬ 
tions with compact support. For any fixed real number x > 0 and 5 > 0, 
set 

{ - for t £ (x,x + S), 

6 (3.20) 

0 for other case. 

In (11.21) first letting gW (t) = <j(£), g^ 2 \t) = 0 and then letting (t) = 0,g^(t) = 
s{t), we have 


1 

6 




'-\-5 

F[k\ ( t , ip)dtdp. 


Since 


and in Z 


i rX+S 

Jim -1 mil = f(x) 


1 r-x+8 

lini u) (p) - J ip [k] (t, ip)dt = u k f (p)<p[ k ] (x, ip), 


we prove the first part of m by the dominated convergence theorem if / £ 
(Co[0, oo)) 2 . For the case of / £ (K 2 (0, oo)) 2 we can approximate / by the 
functions in (Cq[0, oo)) 2 . The second part of (11.31) can be proved similarly. 


4 Proof of Theorem 2 


First let us prove Theorem 2 for a special case. Recall that S and S are the 
solutions corresponding to P = 0 in (II .4|) and (ESI, respectively. 


Lemma 4.1. For f,g £ (L 2 (0, oo)) „ it holds that 

pOO 


-j nOO ^ i f‘ OO 

f(x)g{x)dx = - / Q f (p)Q g {p)dp = - Q f (p)Q@ g {p)dp 

^ J —oo J —OO 


and for x > 0 


f(x) = 


1 


Q f {p)S(x,ip)dp = 


S{x,ip)O f {p)dp, 
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where Qf(p) and Q s (p) are defined by 


®f(p) = / f(x)S(x,ip)dx, @ g (p)= / S(x,ip)g(x) dx. 
Jo Jo 


(4.21) 


Proof. Since it is easy to find that 


S(x,ip) = Q cosh(ipx) + BQ smh(ipx), S(x,ip) = Q cosh(ipx) — QB sinh(ipx), 

(4.22) 

we have 

r°° 1 - i - 

0/(p) = / f(x)S(x, ip)dx = -f{p){Q - BQ) + -f{-p){Q + BQ), 


1 


1 


e s(p) = J S(x,ip)g(x)dx = -(Q + QB)g{p) + -(Q - QB)g(-p), 

where f(p) = / 0 °° f(x) exp(— ipx)dx denotes the Fourier transform of f(x). 
Therefore, by the well-known Parseval equality 

pOO 1 p oo 

/ f(x)g(x)dx = — / f(p)g(-p)dp 


and the identity for u, v £ L 2 { 0, oo) 


u(p)v{p)dp = 0, 


we easily obtain (note that Q 2 = Q) 


^ roo 

7r 


poo 


e /(p) 0 s(p) d P = - / Qf(p)Q®g(p)dp = / f(x)g(x)dx. 

J—oo J 0 

On the other hand, since for all u £ L 2 ( 0, oo) and x > 0 it holds that 

/ oo poo 

u{p) exp(—ipx)dp = / u(—p) exp(ipx)dp = 0, 

-oo J —OO 

we have 

i r°° 

Q f (p)S{x,ip)dp 


^ roo 
7r 


= ^/0*0(Q - BQ)(Q - QB) + ^f(x)(Q + BQ)(Q + QB) = f(x). 

If 00 

Similarly, we can show that — / S(x,ip)Sf(p)dp = f(x). 

,/ — rvi 


□ 


If we put 


pOO 

f{x) = F(x)R{P, 0){x)+ / F(t)I<{P,0;Q)(t,x)dt (4.23) 

«/ X 
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and 


g{x) =R{0,P)(x)G{x)+ f K T (-P T ,0;Q T ) {t,x)G{t)dt, (4.24) 

where F and G can be obtained by solving the above Volterra equations of 
the second kind, then it follows from changing the order of integration and the 
transformation formulae (12.221) and (I2.26[) that 

$/(p) = e F (p), * g (p) = e G (p). (4.25) 

Furthermore, we have 

Lemma 4.2. For f,g £ (P 2 (0, oo)) 4 , it holds that 

p OO pOO p OO p oo 

/ f{x)g{x)dx = / P(x)G(a;)d;r + / / F(y)$(x, y)G(x)dxdy, (4.26) 

J 0 Jo Jo Jo 

where $(x,y) is defined as follows: 

R(P,0)(y)K T (~W,0;W) (x,y) 

+ J ' K(P,0;Q)(y,t)K T (-W^W) (x,t)dt, 0 < y < x, 


$(x,y) = < 


K(P,0;Q)(y,x)R(0,P)(x) _ 

+ J K(P,0;Q)(y,t)K T (-P T ,0;Q T 'j (x,t)dt, 0 < x< y. 

(4.27) 

Proof. On one hand, since R(P, 0)(-) = P _1 (0,P)(-), we have by changing of 
the order of integration 

poo 

I f(x)g(x)dx 

p OO 

F{x)R(P, 0)(®) + j F{t)K{P,0-Q){t,x)dt l 
x |p(0, P)(x)G(x) + J K T (-W, 0;Qr) (t, x)G(t)dt J dx 


fOO /*oo 


J F(x)G(x)dx + J J F(x)R{P,0){x)K T (-P T ,0;Q T ) (t,x)G(t)dtda 

poo pt 

+ / / F(t)K (P, 0; Q )(£, x)R( 0, P)(x)G(x)dxdt 

</o Jo 

poo poo poo --—-- - - 

+J jI J F{t)K{P,0-,Q)(t,x)K T (~-pr,0-,QT) (s,x)G{s)dtdsdx. 
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On the other hand, by (Id.271) . 


pOO /*oo 


nOO 

/ F(y)$(x,y)G(x)dxdy 
Jo 

F(y){K(P,0;Q)(y,x)R(p,P)(x) 


/ 0 JO 

r °° 


/o JO 


K(P,0;Q)(y,t)K T (-P T ,0;Q T ) (a:, f)dt}G(x)dzdy 


fOO poo 


10 Jy 


+ 


F(y){R(P,0)(y)K T (-P T ,0;Q T ) (*,|/) 


A'(P 0; Q){y, t)I< T (-P T , 0; Q T ) (s, t)dt}G(x)dxdy. 


Therefore, to prove (14.261) . it is equivalent to show 


OO poo pOO 


' 0 J x J x 
noo py nx 


F(t)K(P, 0; Q)(t, x)K T (-P T ,0;Q T ) (s,x)G(s)dMsdx 


'0 JO JO 
/>oo /»oo f>y 


F(y)K(P, 0; Q)(y, t)K T (-P^,0;Q r ) (x,t)G(x)dtdxdy 


F(y)K(P,0;Q)(y,t)K T (-P T ,0;Q T ) (x,t)G(x)dtdxdy, 


J0 Jy JO 

which can be easily proved by changing of the order of integration. 


Remark. It follows easily from (12.41) and Lemma 2.1 that 3(v) £ (C* 1 (17)) 

and $(•,■) € ^O^AQ)) 4 . 

Lemma 4.3. For $(x,y) defined by it holds that 


and 

where 


y)B + y) = 0 


$(x,0) = J(x), 3(0 ,y) = C(y), 


J{x) = K T (-P r ,0;Q T ) (®,0), C{y) = K(P,0-Q)(y,0). 
Moreover, the following relation holds: 

J{x) - BJ(x)B = C(x) - BC(x)B. 


(4.28) 

(4.29) 

(4.30) 

(4.31) 
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Proof. For y < x, in view of (12.71) ^ (12.91) in Lemma 2.1, we see by integration 
by parts that 


| ( X ,y)B + B°* {x , y) 

= {R(P,0)(y)Kj (*,!/) 

+ J" K(P, 0; Q)(y, s)K J (-P^ 0; W) (*, s)ds}l3 
+B{R'(P,0)(y)KT (-W, 0;W) {x,y) 

+K(P,0;Q)(y : y)K T (-W, 0;W) (*,»)} 

+B{R(P,0)(y)K£ (-W,0;Q^) (x,y) 

+K y (P, 0; Q){y, s)K * (-W, 0; W) (x, s)ds} 

= {BR'(P,0)(y) - R(P,0)(y)P{y) + BK{P,0-,Q)(y,y) - K{P,0;Q)(y,y)B} 
xK T (-W, 0;W) (x, y ) + K(P, 0; Q)(y, 0)BK x (-W, 0; W) (x, 0) 

+ [ V {. BK y (P , 0; Q)(y, s) + K S (P, 0; Q)(y, a)B - K(P, 0; Q)(y, s)P(s)} 

Jo 

x K T (-W,0;Q^ (x,a)da 

= 0 , 

where we have made use of the relation: B = QB + BQ. For the case x < y, 
the proof of (14.281) is similar. On the other hand, (14.291) is obvious by (14.271) . 

Furthermore, it can be directly verified that the unique solution of problem 
(14.28|) and (14.291) is 

f \{J( X + y)+J(x-y)}~ i B{J(x + y) — J(x — y)}B , y < x, 
5{x,y) = < 

[ \{£( x + y) + c (y - x)} - )^B{C{x + y)~ C{y - x)}B, x < y. 
Consequently, (14.311) follows from the continuity of $(x,y) at x = y. □ 
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Now we apply Lemma 4.3 to show 


Lemma 4.4. It holds that 

®j{p)Q = ®j(p) = ®c(p) = Q®c(p)- 
Proof. By (14.211) . we have 

/•OO f‘00 

®j(p)= / J(x)S(x,ip)dx, Q c (p)= / S{x,ip)C(x)dx, 

Jo Jo 

where S(x,ip) and S(x,ip) are given by (14.221) . Since Q 2 = Q, it is sufficient to 
prove that for all x > 0 

J{x)Q = QC(x), J(x)BQ = -QBC(x). 

First, multiplying right (14.311) by Q, we obtain by QB + BQ = B and C(x)Q = 
C(x) which follows from (12.241) and (14.301) that 

{J(x) - BJ(x)B}Q = C(x)Q - BC(x)(B - QB) = C{x). (4.32) 

Second, since it follows from (12.271) and (14.301) that QJ{x) = J{x), we have 
QBJ(x) = (B — BQ)J(x) = 0. Consequently, it follows from (14.321) that 

QC(x) = Q{J(x) - BJ{x)B}Q = J(x)Q - QBJ(x)BQ = J(x)Q. 

On the other hand, multiplying left (14.321) by B : we have by B 2 = E that 

BJ(x)Q - J{x)BQ = BC(x) = (QB + BQ)C{x), 

that is, 

J{x)BQ + QBC(x) = B{J{x)Q - QC{x)j = 0. 

Thus the proof of Lemma 4.4 is completed. □ 

Proof of Theorem 2. Let C„{x) = "/ a (x)C(x), J„(x) = ') a {x)J{x), where 
the scalar function ^(x) is dehned by (13.31) . It is obvious that both C a and 
J a are continuously differentiable matrix-valued functions with compact sup¬ 
port. Then it follows easily from Lemma 4.4 that Qj a {p) = 0 c a {p)- Hence, 
combining (11.41) . (11.51) . (14.22|) . QJ(-) = J{-), £(-)Q = £(•) and Lemma 4.1, we 
conclude easily that the following matrix-valued function 

If 00 ~ If 00 ~ 

$a(x,y):=- S(y,ip)Q J<7 (p)S(x,ip)dp=- S(y,ip)Q c Ap)S(x,ip)dp 

^ J —OO J —OO 

satisfies the equation 

U X B + BU y = 0, 

and the conditions 


U(x, 0) = Jr,{x), U(0,y) = C a (y) 
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for all x,y > 0. Therefore, if we define 3^(0, 0) = £(0) = J7"(0), then $ a (x,y) = 
$(x,y) in the domain 0 < x,y < a, since 7 a (x) = 1 on [ 0 ,cr] and the two 
matrix-valued functions satisfy the same boundary problem as that in Lemma 
4.3. Moreover, if /, g £ (K^(0,oo)) 4 , then it follows from (14.231) and (14.24[) 
that F(x) = G(x) = 0 for x > a. Consequently, it follows from (14.211) . (14.251) . 
Lemma 4.1 and Lemma 4.2 that 


f(x)g(x)dx 


/•OO /»oo /*00 

/ F(x)G(x)dx + / / F(y)${x,y)G(x)dxdy 

1 0 Jo Jo 

/•OO /*00 pOO 

/ F(x)G(x)dx+ / / F(y)3 a (x,y)G{x)dxdy (4.33) 

/0 Jo Jo 


®f(p){Q + © j a (p)}0 G {p)dp 


®f(p){Q + ®cAp)}®g(p) d P- 


Now define 


D(p) = lim {Q + QjAp)} = lim {Q + &c^(p)} (4.34) 

cr—>-00 cr—>00 

where the limits exist in the sense of convergence of distributions. Indeed, by 
(14.211) and (14.221) we see that both 0 j a ( p ) and Qc a (p) ar e linear combination 
of the Fourier cosine and sine transform of some matrix-valued function with 
compact support. Then it follows from the property of the Fourier transform 
(see e.g. Page 105 in [2T]) that Qj„(p) —> ®j{p) and Oc^(p) —> @c(p) as 
cr —► 00 in the sense of distributions, whence D(p) £ (-Z 7 ) 4 - Therefore, by the 
definition (14.211) we see 

d(p ) = ~{Q + ®j{p)} = ~{Q + e £ (p)}. 

7r 7r 

Thus we can prove the Marchenko-Parseval equality (O) similarly to m- 
Moreover, if one lets g(t) = <;{t)E or f(t) = s(t)E where ?(£) is defined by 
(13.201) . then he can prove (11.81) similarly to (11.31) . □ 
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